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Abstract 

The original purpose of component-based development was to provide tech¬ 
niques to master complex software, through composition, reuse and parametri- 
satiou. However, such systems are rapidly moving towards a level in which 
software becomes prevalently intertwined with (continuous) physical processes. 
A possible way to accommodate the latter in component calculi relies on a suit¬ 
able encoding of continuous behaviour as (yet another) computational effect. 

This paper introduces such an encoding through a monad which, in the com¬ 
positional development of hybrid systems, may play a role similar to the one 
played by 1+, powerset, and distribution monads in the characterisation of par¬ 
tial, nondeterministic and probabilistic components, respectively. This monad 
and its Kleisli category provide a universe in which the effects of continuity over 
(different forms of) composition can be suitably studied. 

Keywords: Monads, components, hybrid systems, control theory 


1. Introduction 

1.1. Motivation and objectives. 

Component-based software development is often explained through a visual 
metaphor: a palette of computational units, and a blank canvas in which they 
are dropped and interconnected by drawing wires abstracting different compo¬ 
sition and synchronisation mechanisms. More and more, however, components 
are not limited to traditional information processing units, but encapsulate some 
form of interaction with physical processes. The resulting systems, referred to as 
hybrid [l, 2}, exhibit a complex dynamics in which computations, coordination, 
and physical processes interact, become mutually constrained, and cooperate to 
achieve specific goals. 

One generic way of looking at components, proposed in Q, emphasises an 
observational semantics, through a signature of observers and methods, that 
makes them amenable to a coalgebraic characterisation as (generalisations of) 
abstract Mealy machines. The resulting calculus is parametric on whatever 
behavioural model underlies a component specification. This captures, for ex¬ 
ample, partial, nondeterministic or probabilistic behaviour of a component’s 
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dynamics by encoding such behavioural effects as strong monads 3 — a per¬ 
vasive mathematical structure with surprising applications in different areas of 
Computer Science (see e.g., 000,10). 

Indeed, each monad captures a specific type of behaviour, which is then 
reflected in the corresponding component calculus. For example, maybe monad 
(1+) introduces partial components; the powerset (IP) monad nondeterministic 
ones; and distribution monad (D) brings (discrete) probabilistic evolution into 
the scene. Can continuous behaviour, prevalent in hybrid systems and control 
theory, be encoded in a similar way, as (yet another) computational effect? Such 
is the question addressed in this paper. 

Monads first came in contact to Computer Science in the 80’s, when E. Moggi 
proposed their use to structure the denotational semantics of programming lan¬ 
guages 0, 0 ■ Later the concept was introduced in programming practice by P. 
Wadler [6(, leading to a rigorous style of combining purely functional programs 
that mimic impure (side-)effects. The key idea is that monads encode in abstract 
terms several kinds of computational effects, such as exceptions, state updat¬ 
ing, nondeterminism or continuations. Such effects are represented by a type 
constructor T (an endofunctor over a suitable category) so that computations 
producing values of type O are regarded as terms of type TO. In this way values 
and computations are explicitly distinguished and programs can be thought of 
as arrows I —> TO representing the computation of values of type O from values 
of type I , while producing some effect described by T. Or, putting it in a differ¬ 
ent way, output values are encapsulated (or embedded) in the effect specified by 
T. A monad comes equipped with an identity and an associative multiplication 
which, from a computational point of view, builds a (trivial) computation from 
a value, and flattens nested effects, respectively. Furthermore, if T is strong 
0 additional machinery is available to distribute the computations’ effect over 
context. The monad structure allows program composition by handling the un¬ 
derlying computational effect through functor T and the flattening operation. 
Actually, each monad gives rise to a so called Kleisli category in which one may 
study the effects of the behavioural type (as specified by the monad) over dif¬ 
ferent forms of composition; ultimately, this leads to rich component calculi (as 
discussed in i)- 

The current paper introduces a (strong) monad K that subsumes the typ¬ 
ical continuous behaviour of dynamical, and hybrid systems. Intuitively, the 
type effect of df (i.e., the underlying endofunctor) represents the (continuous) 
evolution over time of some value in O ; the identity defines a trivial evolution 
(i.e., with duration zero), and the flattening operation allows the control of an 
evolution to be passed along different systems. 

Moreover, the paper explores the corresponding Kleisli category as the math¬ 
ematical space in which the underlying (continuous) behaviour can be isolated 
and its effect over different forms of composition suitably studied. As we will 
see in the sequel, such a category gives rise to several forms of composition 
operators (e.g., sequential, parallel execution), wiring mechanisms, and syn¬ 
chronisation techniques. Again this parallels the role that the categories of 
partial functions, relations and stochastic matrices have as reasoning universes 
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for component composition under the behavioural model provided, respectively, 
by monads 1+, 7 and D [Tl], |T2|- Similarly, this work paves the way to the 
development of a coalgebraic calculus of hybrid components in the spirit of [1]. 


1.2. A tribute to Jose Nuno Oliveira. 

The idea of regarding continuity as a computational effect, or more rigor¬ 
ously, a physical one, entailing a suitable notion of composition and a reasoning 
universe, in the form of a Kleisli category, owes much to the way Jose helped us 
to approach computational phenomena. 

Building on the role of monads in functional pro gram ming and program cal¬ 
culi, as monadic inductive and coinductive schemes |13j, |l4j], Jose introduced us 


to monads both as a powerful structuring mechanism and a source of equally 
powerful genericity. An obsession for patterns and a sharp intuition for generic, 
conceptually reusable structures remain, after all, the hallmark of his illuminat¬ 
ing, socratic teaching. 

In the late 90’s, Jose supervised the PhD work of the second author on 
the coalgebraic calculus of state-based components mentioned above [j]. This 
emerged from the conjunction of two key ideas; first, that a ‘black-box’ char¬ 
acterisation of software components favoured an observational , essentially coal¬ 
gebraic, semantics; second, that the envisaged calculus had to be generic, in 
the sense that it should not depend on a particular notion of component be¬ 
haviour. Monads, actually strong monads, were quickly identified as a source of 
such a genericity, the whole work boiling down to a calculus of monadic Mealy 
machines. Software components were thus studied as coalgebras (in a suitable 
category) typed as 


s->t(sx o y 

where S represents the (internal) state space, and I, O are respectively the input 
and output spaces. T is a strong monad that captures the intended behavioural 
effect. 

Being generic entailed the need for an equally generic reasoning framework. 
By then, the adoption of a pointfree , essentially equational, calculational proof 
style, thus avoiding the somehow more standard coinductive proofs through the 
explicit construction of bisimulations, was understood as the price to be paid 
for genericity, as component laws were to be verified without fixing the working 
monad completely. Generic proofs performed in this style are clear and easy to 
follow, even if often long due to the systematic recording of almost all elementary 
steps. 

For Jose, however, the way proofs are written is not a technicality. Proofs, 
as he taught us every day, are basically honest explanations, bearing evidence 
in a fixed formal context, and therefore must be conveyed in a crisp, clear, 
easily reproducible style, letting the underlying structure to emerge and helping 
to build the correct intuitions. Years later, in the context of a joint research 
project mil ] . Jose championed the use of calculational, pointfree reasoning as a 
way of reinvigorating the role of proof in elementary mathematical education. 
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The pointfree style adopted in many proofs of this paper is also intended as a 
tribute to this view. 

For Jose being generic does not mean to seek refuge in some sort of formal 
ivory tower, of stylised constructions polished ahead of any meaningful intu¬ 
ition. This explains why, being a devoted functional programmer, who resorts 
to Haskell as a pocket calculator, Jose soon started to focus his attention on the 
rich universes of specific monadic computations — their Kleisli categories. If 
pure functions are computations for the identity monad, relations and matrices 
play a similar role in such richer contexts. To be added, of course, and in a very 
concrete way, to the relevant calculator. His systematic, calculational, ‘syntax- 
driven’ work on relation algebra 116. 11711 . as a framework for nondeterministic 
computations, and linear algebra [18|,ll9(, for probabilistic ones, was responsible 
for a fresh understanding of the Kleisli categories of two fundamental monads, 
and lead to a number of new results and simpler, elegant renderings of old ones. 
Having introduced a monad for continuity, this paper initiates the unravelling 
of the corresponding Kleisli category, as the reasoning universe for continuous 
processes, thus, and once again, pursuing a path Jose will certainly cheer. 

1.3. Document structure 

After a brief detour on preliminaries and notation in Section [2] the contin¬ 
uous evolution monad (Jf) is introduced in Section [3] In Section [2 we explore 
the corresponding Kleisli category: as we will see, its arrows define continuous 
systems I —» 3~CO (technically, preliminary versions of dynamical, and hybrid 
systems) and (Kleisli) composition makes possible for a component to execute 
after another, starting its evolution when the preceding one finishes its own. In 
Section [5j we take advantage of the so called Kleisli adjunction to define wiring 
mechanisms and characterise (co)limits. The latter give rise to new forms of 
component composition and corresponding laws. In order to add synchronisa¬ 
tion techniques to our (monadic) framework, Section [G] provides extra structure 
to the underlying functor of monad K. After this we suggest a feedback opera¬ 
tor. In Section[7l we show that monad Jf is strong; this brings us closer to hybrid 
systems as coalgebraic components (in the spirit of I) whose behavioural effect 
is captured by JJ. Formally, coalgebras typed as 

S-^M(SxO)'. 

Finally, Section [8] discusses related work, provides possible research directions, 
and presents concluding remarks. 

In order to illustrate the developments of the ensuing sections, a number of 
classical examples of continuous and hybrid systems will be explored under the 
light of the framework reported in this paper. 
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2. Preliminaries 


2.1. Continuous systems 

Technically, we qualify as continuous a system whose output, for any given 
input, is a (continuous) evolution over time; i.e., an arrow typed as 

i—► U o Td 

dG[0,oo] 


where J, O are, respectively, input and output spaces, 0 Td the space of contin¬ 
uous functions Td —> O (the evolutions ), and Td stands for {r £ R>o | r < d}. 
Actually, this definition includes the family of continuous dynamical systems 
that interpret the non-negative reals (i.e., R>o, here denoted by letter T) as a 


time domain (cf. |2y, 21]). Formally, the latter are characterised as functions, 


$ : A' x T A 
A$ : A' —>■ A T 


such that for any t £ T, x £ A 

<F (x, 0) = x (1) 

$(x,t 1 +t 2 )=®(®(x,t 1 ),t 2 ) (2) 

From a monadic perspective, continuous dynamical systems (in the form A$ : 
A —> A' T ) may be seen as programs whose behavioural effect subsumes some 
form of continuous evolution over time. Indeed, as we will see later in the 
paper, such systems are part of a broader family of arrows that live in the Kleisli 
category of monad TC (Top^). In general, law|T]will be an important part in the 
characterisation of Kleisli composition. We will also see that the traditional view 
of hybrid systems - as a family of dynamical (or continuous) systems indexed by 
a (discrete) state space - coincides with ours; and, moreover, that such systems 
also live in Top^ (due to the machinery that makes Jf strong). 

2.2. Notation 

The key role that continuity takes in this work, suggests the category Top of 
topological spaces and continuous functions as a suitable working environment 
for developing the envisaged results. 

In the sequel, whenever the context is clear, a topological space will be 
denoted by its underlying set. Topological spaces A x Y, X + Y correspond to 
the canonical product and coproduct of A, Y, respectively. Also, for any A C Y, 
assume that A has the subspace topology induced by Y. Finally, whenever Y 
is core-compact (cf. jUJ), space A 5 has the exponential topology. 

Category Top is (co)complete; this allows to take advantage of isomorphisms 
a : (A x 7) x Z = A x (7 x Z), and sw : A x Y = YxX. Top also provides 
a set of useful rules for showing continuity; Figure Q] sums up the ones used in 
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/ : X xY -> Z 


f '■ X Y,g : Y Z 
g.f:X^Z 

f:X->Y 1 ,g:X->Y 2 
(/, g) :X^Y 1 xY 2 


(•) 

(x) 


Ja-A- fY 


(U) 


(A) 


A f :X^Z Y 

f:X 1 ^Y,g:X 2 ^Y 

[f,g]:X 1+ X 2 ^Y 

f '■ X —> Y, img / C B 






( + ) 


4) 


with f a = f • i (for i:A4l) with l • f B = / (for l : B ^ Y) 


Figure 1: Continuity rules in Top. 


the paper. In rule (A), y must be core-compact so that the evaluation function 
ev : X Y x Y —>- X is well defined (c/. 0). 

Universal arrows X —> 1 to the final object in Top are denoted by !, and a 
function constantly yielding a value x by x. Given two functions f,g:X—*Y, 
and a predicate p, we introduce a conditional expression f <\ p g : X —> Y, 
defined by, 


(/ < P > g) x = 



if p x 
otherwise 


Whenever found relevant, and no ambiguities arise, we will denote expression 
(/ Op > g) x by (/ x O p x O gx). The continuous functions minimum 
X : T x [0, oo] —> T and truncated subtraction 0 : T x [0, oo] —> T play a key role 
in some proofs. They are defined by the following equations 


X = TTl O (<) O 7T 2 
© = (-)< (>) > 0 


where <, > are the usual ordering relations over the reals with infinity. 

As usual, functions 7Ti : X x Y —> X, n 2 : X x Y —> Y correspond to the 
projections associated with any binary product, and i± : X —> X + Y, i 2 : 
Y —> X + Y the coprojections associated with any binary coproduct. Moreover, 
symbol * is used to denote the element of a singleton set, and |C| to represent 
the class of objects of a category C. Finally, to avoid a burdened notation, we 
will often drop the subscript in a component of a natural transformation. 


3. The continuous evolution monad 

As mentioned above, we regard continuous systems as arrows of type 

I—> o Td . 

d£[ 0,oo] 
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In order to define them in Top, we need to equip the target object with a 
suitable topology. A first choice would be the coproduct topology (as suggested 
by the expression above), but this is not suitable, since in many cases such 
a topology forbids the system to change the duration of its evolutions along 
different inputs. 

Let us thus explore an alternative topology; the strategy will be similar to 
the one used in the definition of a Moore path category where, given a topo¬ 
logical space X, arrows are paths (i.e., evolutions) [0, d] -X X and composition 
corresponds to the concatenation of those paths (c/. [23j]). Actually, the flat¬ 
tening operation of monad df, discussed below, can be seen as a more general 
version of path concatenation. 

Consider, with no loss of generality, that all evolutions have domain T. Such 
is possible when one notices that T d (for some d £ [0, oo]) is a retract of T 
through the truncation function (the retraction) 

Ad : T —> T d 

X d = id <\ (< d ) > d and considers just those functions / £ O t that become 
constant after time instant d, i.e.f • X d = f. This gives a family of bijections 
{/ £ O t | / • X d = /} = 0 Td indexed by durations d £ [0, oo]. Continuous 
systems thus become arrows typed as, 

i—>{(f 1 d)eO J x[0,oo]\f-x d = f} 

where the target object comes equipped with the canonical topology. This leads 
to the following definition for the underlying functor of monad df. 

Definition 1. df : Top —>- Top is a mapping such that for any objects X,Y £ 
|Top| and any continuous function g : X —> Y, 

dfx = {(f,d)ex j xD\f-x d = f} 

“Mg = g T x id 

where D = [0,oo] is the one-point compactification o/K>o (cf. 0), and g T f = 

9 • /• 

Theorem 1. df is a functor. 

Proof. We need to show that for any continuous functions g : X T, h : Y —> Z, 
Kg : tKX 3~CY is continuous, and df(/i • g) = df/i • 3~Cg 

Since tKg = g T xid , and g is continuous, then “Mg must be as well. Distributivity 
of composition follows from property 

f Kfl = ((- xD).(_) T 3)-^ 

where l is the inclusion map XCX (X T x D), ( _ x D) is the (D) product 
functor, and ( _ ) T the (T) exponential functor. □ 

Let us explore some examples of continuous systems characterised as arrows 
I -x dfO. 
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Example 1. Signal generators are classical examples of continuous systems 
that can generate sinusoidal waves as output. They can be regarded as arrows 
s : M —> !KR such that sr = (r + ( sin _), oo). 

Note that, in contrast to the coproduct topology (in the target object), the 
topology chosen for “K allows durations to change, and thus captures a wider 
range of behaviours. For example, 

Example 2. Consider a thermostat c : R —> J{R that, given the current tem¬ 
perature, linearly raises it to, say, 20 ° C. Such a behaviour can be expressed as 
cr = ((r + - ), 20 © r) where 0 : R x R —> R is the truncated subtraction, i.e., 

e = (-) <1 (>) > o. 

The execution time of system c is thus inversely proportional to the current 
temperature (which is given as input). 

One may also consider another component that takes action after c, and 
whose functionality is, for instance, to maintain the current temperature. The 
result is a composed system that can raise temperatures to a desired level and 
then maintain them - we will explore this specific case in the next section. Of 
course, analogous behaviour can also be found in e.g., cruise control systems, 
water level regulators, and production lines. For example, imagine a component 
of a cruise control system that gives control of the car’s velocity to another 
component whenever an obstacle is detected, or the emergency mode becomes 
active. As we will see in the sequel, Kleisli composition (for monad d£) caters 
for this sort of action. 

The following definition will help in the development of monad “K. 

Definition 2. For any given topological space X € |Top|, define continuous 
function Ox '■ JfX —> X such that 

9 X ( f,d)= f 0. 

Actually, we can canonically extend Ox '■ “KX —> X to a natural transformation 
0 : J~C —y Id, since it is straightforward to show that the following diagram 
commutes for any continuous function / : X Y. 



Moreover, it becomes possible to express the first law of continuous dynamical 
systems (recall the previous section) in a concise, diagramatic manner: simply 
by saying that system c : I —> “KI obeys the first law ((1) above) iff the diagram 
below commutes. 



I 



Actually, we can generalise the diagram to 



I 


where t : / is the inclusion map I' C I. We qualify as pre-dynamical any 

system that follows this generalised condition. Note that both examples above 
dT] and [2]) concern pre-dynamical systems. 

We shall now discuss how to equip Of with the structure of a monad. As 
already mentioned, in programming semantics a monad captures a behavioural 
effect and provides mechanisms to wrap a value into such an effect and to 
flatten two effects into a single one. Technically, they are referred to as the 
monad identity rj : Id — > J~C, and its multiplication /i : 9~CK —> !K, respectively. 
Let us start by defining the unit operation 77 : Id — > df, which will denote trivial 
evolutions. 

Definition 3. Given a space X £ |Top|, function r]x ■ X — » 3~CX is defined by 

r/x x= ( x , 0). 

Intuitively, arrow r)x '■ X —> TCX defines a system whose outputs are always 
trivial evolutions, i.e., with duration zero. For this reason we will refer to r/x 
as copyx, and often omit the subscript. 

Lemma 1. The mapping r) : Id tK is a natural transformation, i.e., for any 
topological space X, rjx ■ X —> TCX is a continuous function, and, moreover, 
the diagram below commutes 


X 


-+Y 


r/x 

tKX 


Kh 


r/Y 

-> “KY 


for any continuous function h : X —»■ Y. 

Proof. To see that rjx is continuous, observe first that yx = (A7Ti, 0). Then, 

7Ti : X x T -> X 


Atti : X -> X 1 


(A) 


(Atti,0) : X —> X T x D 


(x) 

Ur) 


(Atti,0) : X tKX 
It remains to show the naturality of y : Id —> TC. Consider the diagram 
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X h 


->■ hx 

VY 


(x. 0) i. > ( h x , 0) 

h r X id 


where h : X —> Y is an arbitrary continuous function. Property h • x = h x 
entails its commutativity. □ 


It is also simple to see that, for any topological space X £ |Top|, the following 
diagram commutes 

X—^XX 

N ^\ i J x 

X 

(i.e., that rjx is pre-dynamical). Actually, this is one of two laws that charac¬ 
terise 9x as an Eilenberg-Moore IK-algebra [25|, a notion we will visit later in 
the paper. 

The next step is to define multiplication /.i : TCK —> 3~C. We start with an 
(auxiliary) definition of evolution (or path) concatenation. 

Definition 4. Given any elements ( f,d),(g,e ) £ 3~CX, define 

(f, d) * ( 5 , e) = (/ *d 9, d + e) 

where f # d g = f <3 (< d ) > g (_ - d). 

Let us omit the subscript in -h- d . Note that / -h- g is continuous whenever the 
endpoint of / and the startpoint of g coincide. We will show that this condition 
is always met for the case of multiplication. 

Definition 5. Given any topological space X £ |Top|, define 

(f d)= I ^ ‘ ^ ^ ^ tfd^oo 

^ X ’ |^(0 •/, 00) otherwise 

Intuitively, multiplication will serve to concatenate the resulting evolutions of 
two components. 

Lemma 2. The family of mappings p defines a natural tranformation. 

Proof. In appendix. □ 


Lemma 3. For every topological space X £ |Top|, the diagram below commutes 


TVKX TCX 


xce x 

:kx■ 


Ox 


Ox 


-4 X 
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Proof. Consider a pair (f,d) € KKX, where d is finite. Then, 

0-d(f,d) 

= { Definition of /i } 

= { Definition of -H- on point 0 } 

0((0-f,d)) 

= { Definition of "K } 

o • ko (/, d) 

Proof for the case in which d is infinite is achieved via an analogous reasoning 
process. □ 

This property, together with the fact that Ox • rjx = id (discussed above), entail 
that Ox ■ KX — > X is an Eilenberg-Moore df-algebra. In words, an algebra of 
functor K that is compatible with the monadic structure defined above. This 
notion will be rather useful in the sequel. 

Theorem 2. (K,r],p) forms a monad. 

Pi'oof. In appendix. □ 

4. ... and its Kleisli category (Top^) 

If a monad abstracts a computational effect, its Kleisli category, represents 
the universe of computations encapsulated in such an effect. Hence, in the 
case of monad K, the associated Kleisli category of K (Top^) provides an 
interesting setting to study the requirements placed by continuity over different 
forms of composition. Actually, the envisaged calculus of continuous, and hybrid 
components is essentially its calculus. 

This section studies the Kleisli composition of Top^, and illustrates its 
application to the specification of continuous systems - the hybrid ones will be 
discussed later in the paper. We start with the definition of Top M . 

Definition 6. Category Top M is defined as follows: 

• I t °P:kI = I Top |, 

• for any objects 1,0 £ iTop^l, Top ;K (/, O) = Top(/, JfO), and for any 
object I € iTop^l, rji is its identity. 

• Given two arrows C\ : I IKK, C 2 : K KO their composition, denoted 
by C 2 • Ci, is given by p • Kc 2 • C\. Diagrammatically, 
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I 




y'K'KO 


CKc2 



Whenever found suitable, we will denote an arrow c : I —> 3~CO as c : I -+* O, 
and m • c as f c : I —> O t . 

Recall that arrows c : I -+* O are here interpreted as continuous components, 
which means that the Kleisli composition of !K can be seen as a component 
operator. Let us explore its behaviour: consider two systems 

Ci :. I -H- K. C 2 : K —1-> O. 

For a given input x € /, compute the execution time of C 2 • ci, 

7T2 • (C2 • Cl) ( X ) 

= { Kleisli composition } 

7T 2 • H • Jfc 2 • Ci (x) 

= { Definition of let d = 1 x 2 ■ c\ ( x ) } 

7T2 • M(C2 • (/ Cl X),d) 

= { Definition of p } 

d + 7r 2 (c 2 • (f Cl X) d) 

— { Composition } 

d + Tr 2 (c 2 (/ Cl xd)) 

This means that the execution time of C 2 • Ci is the sum of the execution times 
of ci (for input x) and C 2 (which receives value f Cl x d as input). On the other 
hand, 

7Tl • (c 2 • ci) (x) 

= { Kleisli composition } 

7T1 • H • !KC2 • Cl (x) 

= { Definition of !H, let d = • ci ( x ) } 

TTl • d (C2 • (fc 1 X),d) 

= { Definition of fi } 

6 • c 2 • (f Cl x) * (f C2 (f Cl xd )) 

= { Definition of -H- } 

6 • C 2 (f Cl X _ ) < (<d) > fc 2 ( fc 1 Xd)(. -d) 
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Hence, if C 2 is pre-dynamical, 


/(C2 • Cl) if Cl X - ) <i i^d) C> / C2 (/ci x d) (_ d ) 


The last expression tells that for the duration of Ci x, C 2 • ci x evolves first 
according to ci, and then, on its termination, according to C 2 which receives as 
input the endpoint of f Cl x. Clearly, this is the expected behaviour according 
to the definition of operation //, which ‘concatenates’ evolutions. Intuitively, 
C 2 • Ci may also be described as mentioned in Section 1: component ci acts and 
then, at instant d. gives control of its evolution to C 2 . 

If, however, C 2 is not pre-dynamical, then up to completion of interval [0, d], 
C 2 ‘alters’ the evolution of ci; then it proceeds according to its own evolution. 
These notions are illustrated in the following examples. 

Example 3. Given two signal generators ci,C 2 : R -+> R defined as 

Ci r = (r + (sin _), Sir), C 2 r = (r + sin (3 x _), 37 t) 

the evolution c\ • (C 2 • Ci) 0 is represented by the plot below. 


Cl • (C 2 • Cl) 0 



This type of signal is commonly seen in frequency modulation: the varying 
frequency is used to encode information for electromagnetic transmission. Note 
that ci gives control for some time to C 2 , and then ‘takes it back’. 

In order to amplify signals, one can use component a : R -+->• M, where 
a r = ( r x 2 ,0) (note that since system a is not pre-dynamical it can alter 
evolutions of other components). Given input 0, system ci • (C2 • (a • Ci)), 
returns the following evolution. 
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Cl • (c 2 • (a • Ci ))0 



Example 4. Suppose the temperature of a room is to be regulated according to 
the following discipline: starting at 10 °C, seek to reach and maintain 20 °C, but 
in no case surpass 20.5 ° C. To realise such a system, three elementary compo¬ 
nents have to work together: c\ to raise the temperature to 20 ° C, component C 2 
to maintain a given temperature, and component C 3 to ensure the temperature 
never goes over 20.5 0 C. Formally, 

ci x = ((a; + _), 20 0 x ) 

C 2 x = ( x + (sin _), 00 ) 

C 3 x = (x < (x < 20.5) > 20.5 ,0 ) 

In a first try one may compose C 2 , Ci into C 2 • Ci. This results in a component 
able to read the current temperature, raise it to 20 °C, and then keep it stable, 
as exemplified by the plot below. 


C2 • Cl 10 



If, however, temperatures over 20.5 °C occur, composition C 3 • (02 • ci) puts 
the system back into the right track as illustrated in the following plot. 
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C3 • (C2 • Cl) 10 



Clearly, C 3 can be regarded as a supervisor system that, for the sake of effi¬ 
ciency, only acts when temperatures exceed the threshold, using just enough 
power to keep the temperate below the limit. Actually, note that C 3 is able to 
play a supervisory role precisely because it is non pre-dynamical. Of course in 
this specific case, we assume that C 3 has an idealised behaviour, which, despite 
pedagogical, is quite unrealistic. 

The examples above hint at an interesting property of evolutions with infinite 
duration. 

Theorem 3. Consider two arrows C\ : / —H- O, C 2 : O —H- O. If system C 2 is 
pre-dynamical and inrg (7 r 2 • Ci • if C { 00 } for some embedding t : I' I, then 

(C2 • Cl) • 4 = Cl • i 


Proof. 


(C 2 • Cl) • i 

{ Kleisli composition, img (712 • ci • t) C {00} } 

h (C2 • (fa • i)>°°) 

{ Definition of /j, } 

{0 • C 2 • (f ci ‘ if, OO) 

{ System C2 is pre-dynamical } 

(fc 1 • b 00 ) 

{ Notation } 

Ci • i 


□ 

Corollary 1. If c 2 is pre-dynamical and img ( 7 t 2 • ci) C { 00 }, then c 2 • ci = ci. 
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This means that if evolutions of the first component always exhibit an infinite 
duration, the second one, if pre-dynamical, will never have the chance to execute. 

In general, Hf-Kleisli composition provides the basic composition mechanism for 
continuous components; the structure of Top M yields its basic laws. To be more 
concrete, take copy as the trivial system that outputs its input with duration 
zero (i.e., the unit of monad df). Then, given systems Ci,C 2 ,C 3 


copy • ci = ci 

(3) 

ci • copy = ci 

(4) 

(c 3 • C 2 ) • Cl = c 3 • (c 2 • Cl) 

(5) 


5. Wiring mechanisms and (additional) composition operators 

In a category, (co)limits are a main tool to ‘build new arrows from old ones’, 
which in the case of Top^ translates to new forms of component composition. 
Actually, coproducts are easy to obtain through the canonical adjunction be¬ 
tween Top and Top, H . 


L 



R 


which entails that Top M inherits colimits of Top through L. For notational 
simplicity, given a continuous function / : X —> Y, we will denote system 
Lf = V • f ■ X -H> Y by /. 

In Top^, the coproduct (also known as a choice operator) is inherited as follows: 
given two components 


h 



define component [ci, C 2 ] : I\ + I 2 -+* O which makes the following diagram to 
commute. 



Intuitively, [ci, C 2 ] behaves as ci whenever input I\ is chosen, and as C 2 otherwise. 
Such a mechanism is useful to aggregate systems with the same codomain; the 
result being a singular system with different modes of operation (corresponding 
to the respective subcomponents), chosen according to the input received. As 
usual, a functorial sum operator is easily defined. 
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Definition 7. Consider components C\ : I± -+* 0\, c 2 : I 2 -+4 02- Then define 
component C\ EE c 2 : I\ + I 2 -+> 0\ + O 2 as 

Cl EB C2 = [il • Cl, *2 • c 2 ] 


The definition of operator choice as the coproduct universal 

arrow in Top^, 

yields a number of useful laws for free. 


C 3 • [ci,c 2 ] = [c 3 • ci,c 3 • ca] 

(6) 

(ci EB c 2 ) • ii = ii • ci 

(7) 

(ci EB C 2 ) • *2 = *2 • c 2 

(8) 

copyx ffl copy Y = copy x+Y 

(9) 

(d\ EB d 2 ) • (ci ES c 2 ) = (di • ci) S3 (d 2 • c 2 ) 

(10) 

[di, rf 2 ] • (ci EBc 2 ) = [di • ci, g? 2 • ca] 

(11) 

Moreover, 


Lemma 4. For any continuous functions f : X\ —> Y\,g 

: X 2 — y ^2 7 the 

following equation holds 


fBg = f + g 

(12) 


Proof. 

fSg 

= { Definition of E 0 } 

[*1 • f,i 2 • ?] 

= { L is a functor } 

[*1 • f,i2-g] 

= { Definition of L } 

[copy • Zi • /, copy • *2 • 5] 

= { Universal property of coproduct } 

copy • [ii • /, i 2 • y] 

= { Definition of +, definition of L } 

f+'g 


□ 

The left adjoint is also useful to lift functions to the universe of Top K . This pro¬ 
vides a number of interesting operations and wiring mechanisms. For example, 
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recall the diagonal function A: X — > X x X which duplicates the input value; 
the corresponding lifted version A : X -+-»• X X X duplicates evolutions. Take 
now the scalar multiplication * s : ffi. —> R; operation *~ s : K -H> K. can be used to 
amplify signals, a ubiquitous procedure both in signal and control theory. An¬ 
other example is Ff : X x Y -+* X (resp. '■ X x Y -H- X ) which eliminates 
the right (resp. left) side of ‘paired’ evolutions. Finally, sw : X x Y -+->■ Y x X 
swaps the order of evolutions, a functionality graphically represented by wire 
swapping. 

Since L is a functor, the following laws also come for free 

id = copy (13) 

9 • ?=f~f ( 14 ) 

Finding limits in a Kleisli category through left adjoint L is often more diffi¬ 
cult. However, under specific conditions, L also preserves limits. The following 
theorem makes such conditions precise. 

Theorem 4. Consider the Kleisli adjunction L H R of a given monad (T, 77 , y). 
Functor L preserves whatever limits T does. 

Proof. Observe the diagram 



U 


T 

where C 7 is the Eilenberg-Moore category for monad T [ 2 ^] , and K the corre¬ 
sponding ( fully faithful) functor such that T = UKL. Then, consider a limit 
lim^- D in C and assume that T preserves it. This means that T(lim<_ D) is the 
limit of 7D, and equivalently, f7A'L(lim<_ D ) is the limit of UKLD. Since both 
U and K reflect limits, L(lim<_ D) must be the limit of LD. □ 

Note that the theorem above was stated in general terms and is thus applicable 
to any monad. Even though easily proved, its consequences are quite useful. 
For example, in the case of tK it provides pullbacks in Top w , as 

Theorem 5. Functor TC preserves pullbacks. 

Proof. In the appendix. □ 

More concretely, theorems 0] and [5] assert that any cospan A -4 C ■£- B in Top 
gives rise to a pullback in Top^, diagrammatically described as 
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Ax c B —A —> B 
j 


A ■ 


-+C 


! j 

One interesting cospan, worthy of special attention, is A 1 B , which 
induces the pullback 


Ax x B^—>B 
J 


A ■ 


H-> 1 


Indeed, such a construction brings parallelism up front, and moreover, makes 
possible to combine evolutions. More concretely, the diagram states that when¬ 
ever two systems are compatible - in the sense that for any input they produce 
evolutions with equal duration - a new component that encapsulates their par¬ 
allel composition can be defined. Formally, two systems ci : / -+* A, C 2 : I -*->• B 
are called compatible when the diagram 



commutes (note that this is not trivially true, because 1 is not a final object in 
Top^). Then, let E denote set { ((/, d), (g, e)) G KA x KB \ d = e }. When the 
two systems are compatible, a new component {(ci, C 2 }} :/-*->• (A x 1 B) comes 
forward through the mediating arrow (of the pullback), as follows 

((ci,c 2 )) = 7 - (ci,c 2 ) 


where/- ’—U e - >K(Ax 1 B), 7 ((/, d), (g, d)) = (( f,g),d ). 

Note that img (ci, c 2 ) C E precisely because of the assumption of compatibility 
between components (c/. proof of Theorem[5]). In order to keep notation simple, 
we will omit the 1 in the subscript of (d Xi 5). 

We call ((ci, c 2 )) the strict parallel composition of ci and c 2 . Let us illustrate 
its behaviour through a number of examples. 

Example 5. Consider two signal generators, 

ci x — (x+ (sin _), 20 ), c 2 x = ( x + sin (3 x _), 20 ) 

For input 0, system ((ci,c 2 )} exhibits the following behaviour 
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«C1,C 2 )) 0 



0 5 10 15 20 


X 


Consider now component + which adds incoming signals. Then, 

for input 0, the composed system + • ((ci,C 2 )) yields the following signal. 


+ • {{ci,c 2 )) 0 



Since strict parallelism comes from a pullback, the following operator arises in 
a canonical way. 

Definition 8. Consider two continuous systems c\ : I± -+* 0\, C 2 : I 2 -*->• O 2 
such that ci • 7 tT and C 2 • fC. are compatible. Then, define ci B C 2 : h x I 2 -*->• 
Oi x O 2 as 


ClMc 2 = ((ci • 7T1, C 2 • 7rJ}} 

Moreover, the following laws come for free, further contributing to an emerging 
calculus of continuous and hybrid components: in each equation below, assume 
that both its sides are well defined (■ i.e. that the compatibility conditions are 
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respected). Then, we have, 

((ci,c 2 )) • d= ((ci • d,c 2 • d)) (15) 

7fi • (Cl E c 2 ) = Cl • 7 tT (16) 

7tJ • (ci E C 2 ) = C 2 • 7rJ (17) 

((ci,c 2 )) = (d E c 2 ) • A (18) 

copyx E copyy = copyxxY (19) 

(di E d 2 ) • (ci E c 2 ) = (di • Ci) IE {d 2 • c 2 ) ( 20 ) 

(di Ed 2 ) • ((ci,c 2 )) = ((di • Ci, d 2 • c 2 )) ( 21 ) 


Strict parallelism yields a result dual to Lemma 0] 

Lemma 5. For any continuous functions f : X\ —> Y \, g : X 2 —> Y 2 , the 
following equation holds 

f^g = fxg (22) 

Proof. 

f®9 

= { Definition of Kl } 

((/ • ^i,9 • ttJ)) 

= { L is a functor } 

((/ • 7Tl,5^)) 

= { Definition of L, Xi (in Top^) } 

7 • (V f ' ni,V • 9 4 7r 2) 

= { Universal property of product (in Top) } 

7 • (V X 7?) • </ • 7Tl,5 • 7T 2 ) 

= { 7 ' (Wi *Vy 2 ) = Vy 1 xY 2 , definition of X (in Top) } 

V * (/ X ff) 

= { Definition of L } 

/ X 5 

□ 

In some cases, however, putting two components in strict parallel may be too 
restrictive or not enough to meet the system’s design requirements. The next 
section introduces a more relaxed version of parallelism where synchronisation 
comes into play. Mathematically, our construction explores the monoidal nature 
of functor IK. 


21 


6. Synchronised product and feedback 

Synchronised parallelism is a form of composition in which components no longer 
need to be compatible in order to be put in parallel. Instead, each of them can 
change the duration of the corresponding evolutions according to the behaviour 
of the other. The price to be paid is that the previous pullback (or any limit 
in general) is no longer a suitable formalisation. Actually, adding a nronoidal 
structure [26} to functor K, as we will see in the sequel, seems to be a better 
alternative. 

Definition 9. We say that functor K is monoidal (with respect to x) if it 
comes equipped with a morphism m : 1 —> Kl, and a natural transformation 
S : K x K —> K that make the following diagrams to commute for any topological 
spaces X, Y € I Top I • 


(KX x KY) x KZ 

Sxid 


■> “KX x (KY x KZ) 

idxS 


K(X x Y) x KZ 
6 

K((X xY)xZ) 


KX x K(Y x Z) 
s 

■> K(X x (Y x Z)) 


‘KX x 1 


KX <- 


idxm 


J-Ctti 


-4 KX x Kl 
s 

-K(X x 1) 


1 x KX 


KX 


mxid 


0<7T2 


->K1x KX 
s 

- K(1 x X) 


Hence, functor K can be made monoidal once a suitable morphism in : 1 — > Kl 
and a natural transformation <5 : K x K — »• K are defined. 

Definition 10. Let us define such mappings as 

m = copy 

Sx,y ({f,d),(g,e))= {{f,g),dr e) 

where continuous function Y : D x D —> D is defined as Y = < (>) > 7r2. 

As a side note, observe that a possible definition of S resorts to the minimum 
function X (instead of Y) but then the diagrams above would not commute. 
Indeed, for such an alternative to work, m would need to be changed into a 
variant of function copy whose evolutions are always infinite. 

Lemma 6. S is a natural transformation. 
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Proof. We know that function 8 : 3CX x 3CY —t I K(X x Y) is defined as, 

3CX x3dY — )I T xfxDxDA(Ixy) T xDxD ^ 3i{X x F). 

Since T : D x D —> D is continuous, <5 : KX x JiY —> 3C(X x Y) must be 
continuous as well. To show that the naturality property holds, we reason 

K(axb).5((f,d),(g,e)) 

= { Definition of ‘K and 5 } 

((a x b ) • (f,g), dY e) 

= { Universal property of product } 

((a • f,b- g), dY e) 

= { Definition of 5 } 

S ((a-/, d),(b • g, e)) 

= { Definition of !K } 

5-{3iax%b){{f,d),{g,e)) 

□ 


We can now state the expected result. 

Theorem 6. When equipped with natural transformation 8 and morphism m, 
Jt is a monoidal functor. 

Proof. In appendix. □ 

The monoidal structure (Jf, 8, m) defines a specific operator for synchronised 
parallelism, which behaves as follows: given two components with the same 
domain c\ : I —> 3d A. ci : I —> 3dB, define 8 • (ci,C2) : I —> 3iA x 3iB —>■ 
3{(A x B), to be denoted in sequel by (|ci, c 2 [). 

System (jci, C 2 D runs C\ and C 2 in parallel; however, if one finishes earlier than 
the other, it is forced to stall its evolution so that both components end at the 
same time. In other words, the duration of the shorter evolution is increased by 
keeping it constant until the longer evolution terminates. 

Again, this form of parallelism is a lax version of strict parallelism, the cost 
being that many laws that hold before are now lost. Nevertheless, the monoidal 
structure of 3i still makes straightforward to show the following properties. 


f x g • flci,c 2 D = d/ • ci, <7 • c 2 |) (23) 

a • flflci, c 2 D, c 3 D = 0 c i, d c 2 , c 3 DD (24) 

Tiq • do, copyj) = c (25) 

7T2 • d copy, c[) = c (26) 


Moreover, we are able to canonically define a new operator, following a path 
similar to the one used to define ES and Kl. 
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Definition 11. Given systems ci : I\ -H> Oi, C2 : In -H> O2, component ci0C2 '■ 
I\ x I<2 -+* 0 \ x O 2 is defined by 

Cl 0 c 2 = dci • 7tT , C 2 • 7rJ D 


The following laws arise from routine calculations 


svj • (c 2 m ci) = (ci m c 2 ) • sw 

(27) 

a • ((ci m C 2 ) 0 C 3 ) = (ci 0 (c 2 0 c 3 )) • a 

(28) 

copyx 0 copyy = copyxxY 

(29) 

f®9 = f x g 

(30) 


Note that strict and synchronised parallel composition behave identically but 
with one exception: in any given execution, the latter increases the execution 
time of a system that finishes earlier than the other. Hence, for compatible 
components both operators behave exactly in the same way, and, therefore, the 
former inherits all laws derived in this section for the latter. 

Next, we introduce iteration for continuous systems. This facilitates com¬ 
ponent specification and, moreover, can be used to express (or detect) Zeno 
behaviour Q. 

Definition 12. Given a component c : X -*->■ X, component c n : X —H- X is 
defined by the ( Kleisli ) composition of c with itself n times. Formally, 

c° = copy , c n = c" -1 • c 

It is straightforward to check that the following equations hold. 

copy n = copy 
c 1 = c 

(jiyn X m 

C n • c m = c n+m 
(c ffl d) n = c n md n 
(c Kl d) n =c n Md n 

Infinite iteration leads to the familiar notion of feedback. 

Definition 13. Let (X, d) be a complete metric space, and c : X -+* X a pre- 
dynamical system; denote the series (7^ "C^a^igN by (sj)igN, and the sequence 
(tti • c l (a;))ig N by (/ l ) ieN . 

Then, assume that for any x € X whenever the series (si)»eN converges 
the sequence is Cauchy. More concretely, its elements get progressively 

closer to each other with respect to the metric, 

d*{g,h) = sup d{g[t), h(t)). 
te t 


(31) 

(32) 

(33) 

(34) 

(35) 

(36) 
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The interested reader will find in f¥\ l more details about this metric. 

c c c 

Finally, define infinite iteration (X -+* X -*->■ X -+* ...) as vc : X -+> X 
where 


7 r 2 • vc (x) = 


oo 

lim^oo 


Si 


if the series (Sj)igpj diverges 
otherwise 


(7Tl • vc (x)) t = 


fk t 

( i yoo f'i ) t 


if t < (772 • vc (x)) 
otherwise 


for k the smallest value such that t < Sk- 

Intuitively, to compute the value at a certain instant (t) in the evolution (771 • 
vc (x)), we need to compose c with itself the necessary number of times for the 
composite ‘to reach that instant’; only then it is possible to extract the value. 
To be concrete, if each iteration of c has two seconds of duration, to calculate 
the value at five seconds in the evolution (771 • vc{x )), we consider the composite 
c 3 and compute the expression (771 • c 3 (x)) 5. 

Observe that, since c is pre-dynamical, the calculated value is not changed 
by additional iterations, i.e. 

( 77 ! • c k (x)) t = (77! • (c • c k ) (x)) t. 


Actually, in the definition above one may forget the assumption of c being pre- 
dynamical as long as it is ensured that the sequence (/i)ieN is always Cauchy. 

The following section gives concrete examples of parallel operators and (in¬ 
finite) iteration at work. The role of feedback in handling Zeno behaviour is 
illustrated as well. 


7. From continuous to hybrid systems 

Having characterised a calculus of continuous components based on the struc¬ 
ture of the Kleisli category of monad TC, the next step is to broaden the picture 
in order to handle systems that exhibit continuous and discrete behaviour inter¬ 
twined. Such is the purpose of this section. A number of examples will illustrate 
the approach proposed here as well as some of the operators introduced in the 
previous sections. 

Our aim is to equip continuous systems with an (internal) state space that 
behaves in a discrete manner. Therefore, arrows become typed as 

Sxl —> S x TCO. 

Intuitively, given a state (s £ S) and an input (i £ I), the component transits 
(internally) into another state and presents continuous evolutions that can be 
directly observed. This gets us closer to the notion of hybrid system, as a 


25 


family of continuous systems indexed by a state space. On the other hand, this 
approach is aligned with the notion of components as coalgebras (as described 
in (111). Actually, our aim is to characterise hybrid systems as coalgebras with a 
discrete (internal) behaviour, and (external) continuous evolutions. 

The cornerstone of this move from continuous to hybrid components is the 
notion of tensorial strength for monad ! H: a natural transformation r : IdxH —> 
H(ldxld) that commutes with the monad operations and with specific monoidal 
structure of the base category (see the formal definition in S)- Indeed, tensorial 
strength allows us to transport such systems to Top M , via composition: 

c S x I —y S x HO 

T-c-.SxI-> H(S x O) 

Definition 14. Given topological spaces X , Y £ |Top| a (right) tensorial strength 
of monad H is the function tx,y '■ X x HY —> H(X x Y) defined by 

tx,y (x, (/, d)) = ((x,f),d). 

Interestingly, function r corresponds to the uniform characterisation of tensorial 
strength for monads over Set (cf. [H|). This entails that all diagrams that 
need to commute do commute, and therefore we just need to show that r is 
continuous. For this, observe that r can alternatively be defined as (\T a ,T b ) : 
X x HY -> H(X x Y) where, 

Ta (Or, (/, d)),t) = (: x,ft) 
r b (:r, (/, d)) = d 

Since T a ,r b are continuous, so is r. 

Corollary 2. Natural transformation r : Id x H —> H(Id x Id) defines a 
tensorial strength for monad H. 

Note that one can also define a natural transformation r; : Hx Id H(Idx Id) 
(known as left tensorial strength for H), via the equation 77 = (Hsw) • r • sw. 
Moreover, a monad is commutative, if the equation below holds. 

t • n =n • r 

This is not, however, the case for monad H. as the following counter-example 
reports. 

Example 6. Recall the two signal generators, introduced in Example 0 

ci x = ( x + (sin _), 20 ), ci x = ( x + sin (3 x _ ), 20 ) 

The application of left, and right tensorial strength to the composed function 
(c\, cf) '■ R —> HCR x yields the behaviours depicted below. 
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T • Tj • (ci, C2) 0 



0 10 20 30 40 


x 



x 


Clearly, r • r; 7^ 77 • r; but note that the plots illustrate an interesting aspect: 
specification r • 77 • (ci, C2) reads ‘first let the component in the left to act, then 
the one in the right’; and conversely for 77 • r • (01,02). Moreover, note that 
each component ‘waits’ for the other by stalling the corresponding evolution. 
This introduces yet another synchronisation mechanism. 

Equipped with tensorial strength r, we may now explore two classical examples 
of hybrid systems from a component-based perspective. We start with the 
bouncing ball system. 

Example 7 . Consider a bouncing ball dropped at some positive height and with 
no initial velocity. Due to the gravitational effect, it will fall into the ground but 
then bounce back up, losing, of course, part of its kinetic energy in the process. 

From this description, one may regard the bouncing ball as a hybrid component 
whose (continuous) observable behaviour is the evolution of its spacial position, 
whereas the internal memory records velocity, updated at each bounce. To 
define such a component we resort to Newton’s equations of motion. 

pos a (v,p, t) = p + vt— ^at , vela {v, t) = v — at 
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from which we can derive the function that, given a positive height and a current 
velocity, returns the time needed to reach the ground; formally, 


, , W2ap+v 2 +v 

ZpOS a {V,p) = --- 

Let us then define the discrete behaviour of the bouncing ball bd : V x P —>■ V 
bd (v,p) = vel g {v,zpos g (v,p )) x -0.5 

where 0.5 is the dampening coefficient. For the continuous part b c : V x P —» ‘KP 

b c = ( pos g ,zpos g ) 

where g = 9.8 (Earth’s gravity). The resulting system is a ball bouncing on 
planet Earth, denoted by b and formally defined as b = t ■ (bd, b c ). Assume that 
the initial state of b is 0. Then, through the iteration operator, and assuming 
five as the initial position one gets, for instance, the following behaviour. 



X 


Analogously, we can define a ball bouncing in the Moon (here denoted by let¬ 
ter c), and compare the behaviour of both bouncing balls by putting them in 
parallel, with the same initial state 0. 



X 
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Note that (fe 3 m c 3 ) 7 ^ (fe[Hc) 3 . An interesting question to pose is about the 
durations that components vb and vc output. Indeed, the intuition is that du¬ 
rations are always infinite (since feedback involves infinite sums), however, due 
to the Zeno effect, the durations that concern this example are actually finite: 
they correspond to the time at which the ball stops moving. Such durations are 
given precisely by the computation of tt 2 • {vb) and n 2 • {vc) with respect to a 
given input. 

Example 8. Alternating pumping systems are often used to regulate the water 
level of reservoirs. Consider one that fills two tanks alternatively in cycles of 
ten seconds, which means that some sort of internal memory is required {to 
remember which was the last tank served). 

Thus, the discrete part Wd ■ S x LS is defined as 

Wd = flip • 7Tl 

where S = {T, _L} is the discrete state space and flip the function that switches 
between the elements. Let us assume that the initial state is T. Then, we define 
the continuous behaviour w c : S x L —7 3~CL 

w c {s,(h,l 2 )) = (/ s (Zi,Z 2 ),10) 

where / T (Zi, h) = {{h + _), h) and f±{h, I 2 ) = {h, {k+ -))■ As expected, the 
pumping system is given by equation w = r • ( uid,w c ), which, for input (0,0), 
yields the following plot. 



X 


On a different note, it is natural to consider that the pump takes some time 
to switch from one tank to the other: for illustration purposes let us assume 
that time to be ten seconds. To simulate such a delay we can define a variant 
of copy, denoted by copy 10 , that always outputs evolutions with duration ten. 
Then, again for input (0,0), system {copyio • w ) 3 outputs 
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(copy 10 • to ) 3 (0,0) 



It is also important to analyse situations in which water flows out. Thus, con¬ 
sider a hybrid system 2 : 1 x 1 —> M(1 x L) (with trivial state space 1) whose 
continuous part 


— ((/ 25 / 2 ), 10 ) 

dictates the rate of water flowing out in each tank, here represented by a clock 
that runs at half the normal speed. Then, we specify the result of w and 
2 acting together in the same set of variables. For this, we define function 
h : (S x L) x (1 x L) —> (S x L) x (1 x 1) where 

h ((s, li, I 2 ), (*,£,?/)) = ((Ml Ox,hOy), (*,*)) 

Intuitively, function h subtracts water in accordance with the rate specified by 
component 2 . For input (0, 0), system h • (w Kl 2 ) yields the plot below. 


(h . (to Hz)) 3 ( 0 , 0 ) 



8. Conclusions and future work 

It is well known that software systems are becoming prevalently intertwined 
with (continuous) physical processes. Such an architecture, however, renders 
their rigorous design (and analysis) a difficult challenge that calls for a wide, 
uniform framework combining the continuous and discrete sides of Mathematics. 
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As a first step towards a component-based framework for hybrid systems, in 
the spirit of Q, this paper showed how continuous evolutions can be encoded 
in the form of a strong (topological) monad. As discussed in Section 1, to 
capture specific behavioural models through monads has been a successful path 
in Computer Science: such was the case of nondeterministic behaviour, and the 
(discrete) probabilistic one; but occurrences in the continuous domain also exist. 
A prime example is the Giry monad [29}, which captures stochast ic p rocesses 
and has been object of study in a number of papers (e.g., [ 1 , 30[ 31 > 32[). Along 
similar lines, monad “K provides a categorial universe for continuous, and hybrid 
systems, where the effects of continuity over different forms of composition can 
be isolated and suitably studied. 

This universe, i. e.the Kleisli category Top M , offers different forms of system 
composition, wiring mechanisms, and synchronisation techniques. For example, 
Kleisli composition lets the control of an evolution to be transferred from one 
system to the other, but also allows evolutions to be dynamically modified 
(as observed in the case of signal amplification). Such behavioural patterns, 
as discussed in Section (3j are often found in systems like thermostats, cruise 
control systems, and signal generators. But more generally, in control loop 
systems - traditionally comprised of a network of digital controllers that manage 
a physical process over time through a feedback loop architecture. In this case, 
the controllers, possessing different functionalities, periodically pass control of 
the physical process among themselves. 

The underlying categorial framework hinted at several composition operators 
(through corresponding universal constructions), and facilitated the elicitation 
of several compositional laws. Throughout the paper, the results achieved were 
illustrated with classic examples of hybrid systems, namely a thermostat, a 
bouncing ball, and a water tank system. 


8.1. Related work 

Hybrid automata 33[ are the de facto formalism for the specification of hy¬ 
brid systems. Roughly speaking, they are a variant of classic automata that 
allows variables to continuously evolve while in a state. This defines the contin¬ 
uous behaviour of an hybrid system, which is then paired with discrete actions 
given by the usual state transitions. Parallel composition of hybrid automata 
proceeds similarly to the classic case, where common labels act as synchronising 
events. Interestingly, in [34j] Bornot and Sifakis introduced additional synchro¬ 
nisation mechanisms that make one system wait for the evolution of the other 
to end, or, on the contrary, force it to finish earlier. This seems to be intimately 
related to whatever monoidal structure is given to functor H. 

During the last years there were also developments concerning the addition 
of new dimensions to hybrid automata: for example, [35| shows how to take 
reaction times into consideration in a compositional setting. In our case, we 
took advantage of dawdler components, like copyio , to introduce such delays. 

The ‘rationale’ underlying hybrid automata is powerful, and highly intuitive, 
but in some cases lacks expressive power: for example, those systems in which 
evolutions can be dynamically changed by some of the components are very 


31 








hard to specify. Moreover, aside from parallel composition, the authors have no 
knowledge of deep developments that concern new compositional operators for 
hybrid automata. 

The industrial tool SimulinlI]], on the other hand, offers a highly expres¬ 
sive component-based language, and is thus closely related to the framework 
proposed in this paper. Indeed, SiMULlNK supports a rich palette of compo¬ 
sitional operators, and computational units. It possesses behavioural patterns 
that involve dynamical alteration of evolutions, delays, and synchronisation. 
Moreover, the transfer of the control of some evolution is not hard to define. 
All this renders SiMULlNK a very interesting tool. The cost is the lack of a clear 
semantics, which impairs formal analysis and the elicitation of compositional 
laws - actually, some recent efforts have been made towards the formal verifi¬ 
cation of SiMULlNK models in alternative tools (c/. [36, 371). In addition, the 


components available are rather limited in what concerns the characterisation 
of their internal memory and respective transition dynamics. 

It would be interesting to study the embedding of (a subset of) Simulink’s 
language into Top^. In principle, Top^ could act as a tool complement, pro¬ 
viding a basis for the formal analysis of (critical fragments) of hybrid systems. 
We stress, however, that we do not aim at emulating SiMULlNK, but rather at 
a suitable coalgebraic framework for hybrid components, where we consider the 
discrete transitions to be internal behaviour, and the continuous evolutions the 
observable part. From this point of view, SiMULlNK is very distant from such a 
line of work. 

There is also a close relation between the work here reported and P. Hofner’s 
algebra of hybrid systems [38]: the latter’s main operator is used to concatenate 
evolutions. Moreover, the algebra possesses secondary operators, like parallelism 
and synchronisation, that are equally available in Top^. Our approach, how¬ 
ever, and differently from P. Hofner’s calculus, is structured around a monad 
that encodes the notion of continuous evolution; this brings up a number of 
canonical constructions and smooths the integration with other behavioural ef¬ 
fects, such as nondeterminism or probabilistic behaviour. 

Finally, a few categorial models for hybrid systems have been proposed along 
the last two decades. For example, document [3f| introduces an institution 
in essence, a categorial rendering of a logic - for hybrid systems, and provides 
basic forms of composition such as free aggregation ( i.e ., parallelism without 
interaction) and interconnection where some attributes and events are shared 
between two systems. Around the same time, Jacobs [ 13 ] suggested an object 
oriented coalgebraic framework where hybrid systems are regarded as coalgebras 
equipped with a monoid action: coalgebras define the discrete transitions, and 
monoid actions the continuous evolutions. Some years later Haghverdi et. al 


41] explored the connection between a formalisation of hybrid systems (close 


to hybrid automata) and open maps. The objective was to provide appropriate 
notions of bisinnilation both for dynamical, and hybrid systems. Composition 


1 http://www.mathworks.com/products/simulink 
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mechanisms, however, were not studied in this context. 


8.2. Future work. 

Our next step is the development of a calculus of hybrid components (as in §|) 
based on monad 3~C and its Kleisli category. The calculus from Q, in its coalge- 
braic spirit, is bisinmlation-based, with bisimulation given as the usual span of 
simulations [28|. The framework that this paper sets, however, offers a promis¬ 
ing basis to explore alternative notions of (bi)simulation for continuous and 
hybrid systems. This has points of contact with the work of Haghverdi et. al in 
[4l|; but note that we use coalgebraic machinery, and follow a component-based 
perspective, which makes possible to study the relation between (bi)simulation 
and (the different) compositional operators. 

A second line of research concerns the development of a taxonomy of con¬ 
tinuous, and hybrid systems living in Top^-. Indeed, as Stauner showed at the 
beginning of the century in his PhD thesis |42j, topologies are useful to elicit a 
number of important properties. For example, the notion of robustness (preva¬ 
lent in control theory) becomes simple to formulate: intuitively, a system is 
robust if small changes in the input lead to very similar evolutions. In Top M , 
since each system has a topological semantic base, one can express how robust 
it is by varying the topology in its source object. At one limit, if the topology 
is discrete, the system is seen as chaotic. At the other end, i.e., if the topology 
is indiscrete, the system must always output the same evolution. 

Actually, the compositional nature that underlies Top M allows us to rea¬ 
son about the robustness of the system at hands through the analysis of (the 
robustness of) its simpler constituents. One disadvantage of this approach is 
that composition in Top^ is strict , in the sense that components with different 
topologies in the connecting points cannot be composed. For example, it is hard 
to put a chaotic component after a robust one. Part of our current research tries 
to relax this condition while maintaining stability, whenever possible. 
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Appendix 

Proof of Lemma [2] The proof is divided in two parts: the first establishes 
continuity of the mappings, the second concerns naturality. Consider the map¬ 
ping fix '■ KKX — > KX; we are going to show its continuity. First we observe 
that fix can be alternatively defined as (Aa, b ) where 

ixid cone 

a = KKX x T-> (X TxT x D) x T-> X 

i = ‘K'KX - 1 Xld - > (X T ) T x D ——■» X TxT x D 

for cone ((f,d),t) = f (t X d,t 0 d). The definitions clearly show that a is 
continuous. For function b we have 

n^xid c 

b = KKX -■> D T x D -f D 

where c (/, d) = (/ d) + d <\ (d ^ oo) > oo. Since the canonical restriction 
(+) • (ev,TT 2 ) : D T x T —>■ D of c is continuous we just need to show that the 
latter is continuous at infinity. Actually, this comes for free once proved that 
given any neighbourhood N D (x, oo] in D of oo we can find a neighbourhood 
V in D t x D of (/, oo) such that c ( V ) C N. 

Consider neighbourhood D T x (a;, oo]. It is clear that c (D T x (a;, oo]) C 
(x, oo] C N. 

Next we show that fi is natural, i.e., that for any continuous function h : X —>• Y 
the diagram below commutes. 


KKX KKY 


MX 

KX 


:Uh 
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First we assume that (/, d ) £ KKX has finite duration, 

H • KKh (/, d) 

= { Definition of l K, fi } 

(9 • Kh ■ f, d ) *{Kh- f d) 

= { Q is natural } 

(h • 9 • f,d) -h- (“Kh • f d) 

= { Definition of IK, composition } 

Vih (9 • /, d) -H- ‘Kh (/ d) 

= { (-H-) is natural } 

Kh((9-f,d)*(f d)) 

= { Definition of Jf, fi } 

Kh • ii (/, d) 

The proof for the case in which (/, d) £ KKX has infinite duration is analogous 
to the above. 

□ 

Proof of Theorem [2] We have to show that the following diagrams commute. 

^K 2 

- >K 

Note that the proof below becomes much more simpler if the evolutions involved 
have infinite duration. 

Let us start with the left triangle. 

M • V (/> d) 

= { Definition of r] } 

M ( (A d) , 0) 

= { Definition of } 

(fl •(/,<*),()) »((/,<*) Q) 

= { Definition of constant } 

(g.(/,d),0)*(/,d) 

= { Definition of -H-, definition of constant } 

CM) 
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For the right triangle we have, 


A* • M if, d) 

= { Definition of IK } 

M (v • f, d) 

= { Definition of 11 } 

(0-V d) 

= { Definition of 77 } 

= { Definition of -H- } 

(8 • V • f, d) 

= { Eilenberg-Moore } 

CM) 

It remains to show that the square commutes. Before giving the formal proof, 
we present the corresponding intuition from a geometric perspective. 

Let us then start by observing that an element in KKX , may be intuitively 
seen as a square, where each column is a function in KX. Then, note that 
multiplication (p : J-CKX —► CKX) keeps just the first row and last column of 
the square, as illustrated below. 


-—* y 

As expected, the intuitive picture of an element in CK 3 X is a cube, 


z 



such that a projection on the z-axis yields an element of J-CKX (geometrically, 
a square as described above). 

Let us now observe that, resorting to multiplication, we can reduce the cube 
into a square. Actually, we can do this in two different ways: via fi : K 3 X —> 
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JCHX, or Jf/i : CK 3 X —> HOOT. In the former case, only the front and right 
surfaces are kept (picture below in the left). In contrast, function J{/i applies 
[i to each projection on the rr-axis, and thus only the bottom and back surfaces 
are kept (picture below in the right). 



Finally, applying p, : 0~CKX —> !KA' to the resulting squares yields the same 
result, 
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More formally, we reason 

M • Jf/x (/, d) 

= { Definition of IK } 

M (M • f, d) 

= { Definition of fi } 

{0 • V • f, d) * (m • / d) 

= { Eilenberg-Moore } 

= { Let f d = definition of -h- } 

df) 

= { 0 is natural } 

{(e-o-f,d)*(e-f , ,d'))*(f , d') 

= { Notation (/ d), definition of -h- , definition of IK } 

= { Definition of /i } 

(JC0.M(/,d))*(/' <0 

= { Definition of /i } 

( • n (/, d)) * (7Ti • fj, (/, d) tt 2 • /i (/, d)) 

= { Definition of IK } 

(0 ■ 7T1 • M (/, d), 7T2 • /r (/, d)) -H- (tti • /x (/, d) tt 2 • n (/, d)) 

= { Definition of [i } 

M 0 (/, d)) 

= { Composition } 

M • M (/> d) 


Proof of Theorem [5] Consider the following pullback in Top 


dxcB 
J 


->■5 


/ 


4C 


where / and g are arbitrary continuous functions. We need to show that 


□ 
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Jt(A xcB) ^4 'KB 


i 

5L4 ■ 


Jif 


Jig 

AKC 


also forms a pullback in Top. 

For this, observe that functor 3~C comes from the composition of functors 
( _ ) T , and ( _ x D), both of which preserve pullbacks. Indeed, they give rise to 
the commuting diagram 



where 7 ((ei, d), (e 2 , d)) = ((ei, e 2 ), d). Let us denote 7 • (ci, c 2 ) by ((ci, c 2 )). 

Since functor df forces specific conditions on evolutions (recall that (e, d) £ 
ACX implies e • Ad = e) some work remains to be done. In fact, we need to show 
that img((ci, c 2 }} C !K(AxcB) whenever img C\ C J-CA, img c 2 C J{B, and Ci, c 2 
make the outer square to commute. In other words, we need to show that, under 
these conditions, ((ci,c 2 )) factors through 1 : 5f(A xq B) (A Xc B) T x D; 
diagrammatically, 

X (<Cl^C 2 )>(^ XcB ^T xD 

L 

J{(A x c B ) 


Consider an element x £ X, and denote ((ci,c 2 )} x by ((ei,e 2 ),d). Since by 
assumption e\ • Ad = ei, e 2 • Ad = e 2 , it is clear that (ei,e 2 ) • Ad = (ei,e 2 ) and 
therefore ((ei,e 2 ),d) G J£(A Xc B). 

□ 

Proof of Theorem [6l We need to show that the following diagrams commute. 


(JfX x dfT) x dCZ 
Sxid 


* Jf A x (IKY x J£Z) 

idxS 


X(X x Y) x 'KZ 
s 

mx xY) x Z) -—-- 


•KX x 1K(Y x Z) 
5 

jf(x x (y x z)) 
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:MX x 1-——> MX x Ml 

*1 S 

KX < - K(X x 1) 

0~ClT i v 7 

We start with the upper square. 


1 x MX -—> Ml x MX 

7l"2 S 

KX < -—- K( 1 x X) 


Ma • 5 • (5 x id) (((ei, di), (e 2 , d 2 )), (e 3 , d 3 )) 

= { Definition of S and ‘K } 

( a • ((ei, e 2 ), e 3 ), ((di Y d 2 ) Y d 3 )) 

= { Definition of product, Y is associative } 

((ei, (e 2 , e 3 )), (di Y (d 2 Y d 3 ))) 

= { Definition of 5 } 

<5 ((ei, di), ((e 2 , e 3 ), d 2 Y d 3 )) 

= { Definition of id x 5 } 

6 • (idx 5) ((ei, di), ((e 2 , d 2 ), (e 3 , d 3 ))) 

= { Definition of a } 

S x (id x (5) • a (((ei, di), (e 2 , d 2 )), (e 3 , d 3 )) 

Then, for the diagram above in the left we reason, and proceed similarly with 
the one in the right. 

Mni • S • (id x m) ((/, d), ★) 

= { Definition of in. 5, and r K } 

(tti • </,*), d Y 0) 

= { Cancellation x, 0 is the identity element (for Y) } 

(fid) 

— { Definition of 7ri } 

TTl ((/> d),*) 


□ 


43 



